Hidden order revealed in quantum oscillations in cuprate superconductors 
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We follow the line of reasoning that hidden broken symmetries are the root of quantum oscillations observed 
in underdoped superconductors and examine the role of bilayer splitting and incommensuration. This is a view 
that eschews the notion of a featureless Mott liquid as the source of complexity. Instead, our view is grounded 
in a conventional Fermi surface and quasiparticles. We show that bilayer splitting and/or incommensurate d- 
density wave order can lead to many interesting results, in particular a splitting of the main frequency of the 
quantum oscillations. 



Introduction: Over the past twenty years of the discovery 
of the high temperature superconductors, the imperfections of 
the materials have hidden many important facts and given rise 
to a set of dogmas that are all pervasive. Inspired by the re- 
cent quantum oscillation experiments in underdoped cuprates, 
we depart from the accepted wisdom and place prominence 
to the notion of broken symmetries [1]. Broken symmetries 
associated with phases of matter dictate the elementary exci- 
tations that determine the macroscopic properties, which offer 
us clear clues about the gross features of a physical problem 
and transcends whether or not the system is strongly corre- 
lated. It is this view that inspired us to propose the picture of 
the d-density wave (DDW) that can potentially unify a wide 
range of disparate features of the cuprates 0. 

The quantum oscillation experiments J5] @ [5] |6] |7) indi- 
cate that in highly underdoped cuprates: (1) the de Haas-van 
Alphen (dHvA) effect, the Shubnikov-de Haas (SdH) effect, 
and the oscillations of the Hall coefficient, all have the same 
underlying cause, the conventional quantization of the Landau 
levels, meaning oscillations that are periodic in the variable 
1 /B, where B is the magnetic induction; (2) the Hall coeffi- 
cient in high magnetic fields and low temperatures is negative; 
(3) given the conventional nature of the oscillation spectra, it 
is difficult to avoid the Luttinger theorem [8, 9| regarding the 
volume of the Fermi surface and the number of charge car- 
riers; (4) there are strong indications that Fermi surface re- 
construction due to a broken symmetry plays an important 
role; (5) there are also tantalizing hints that order may even 
incommensurate ifTOll ; if there is incommensurate DDW or- 
der, this will strongly enhance the tendency to inocmmensu- 
rate charge density wave order (CDW) with twice the ordering 
vector [11]. Thus, indirectly through the intermediate CDW, 
there may be a link between the incommensurabilities of the 
spin density wave (SDW) and DDW; (6) dHvA in the mixed 
state can reflect the properties of the normal state. In a recent 
paper, one of us has addressed some of these aspects from the 
perspective of DDW ifTH . 

Here we shall concentrate on the following: (1) the effect 
of bilayer splitting in YJ^^Gu^Oq+s (YBCO), an unavoid- 
able structural constraint that one must address; (2) possible 
incommensuration of the order parameter J2][l3]; (3) the ro- 
bustness of dHvA in the mixed state reflecting the putative 
normal state; (4) some important unanswered questions. Our 
main analysis is based on a Hartree-Fock theory of the order 



parameter and the Gorkov equations for the mixed state. It 
is the beauty of an order parameter theory that those features 
that are protected by symmetries can be explored in the weak 
interaction limit, although the system may be strongly inter- 
acting. 

Bilayer splitting Bilayer coupling, ij_(k), in YBCO is 
well-known. It has been parametrized in terms of a mo- 
mentum conserving tunneling matrix element. For tetragonal 
structure it is lfT4l[T5l 

i±(k) = -i [cos(fc x a) - cos(k y a)} 2 , (1) 

where a is the lattice spacing. 

The outline of the calculation is as follows. The total 
Hamiltonian H = H i + H2 + H12 , where Hi and H2 are the 
full Hamiltonians of the layers (1) and (2) and their spectra 
are degenerate, as they are identical. The tunneling Hamilto- 
nian H12 is given in the momentum conserving case to be 

tf 12 =^(k)(4 ( >g + 1^2) (2) 

k,<r 

Commensurate order: We now argue that in the Hartree- 
Fock approximation, Hq = H\ + H2 can be written as the 
effective DDW Hamiltonians 

tt _ ( Ml) m f(i) (1) ,i^9\ 

kEB.BZ,a 

+ E (^k c kl 1) 4+ Q , CT + l-2 + h.c.) (3) 

kEB.BZ,a 

where Q = (tt / a,n / a), and is the single particle spectra. 
The superscripts on the electron creation and annihilation op- 
erators stand for the layer index. The reduced Brillouin zone 
(RBZ) is bounded by k y ±k x — ±ir / a. We shall parametrize 

e k by ma 

ek = — 2i(cos k x a + cos k y a) + At' cos k x a cos k y a 
- 2t"(cos2k x a + cos2k y a) (4) 

and the DDW gap by Wk = ^(cos k x a — cos k y a). With 
the choice of the quadratic Hamiltonian Ho in Eq. Q, it can 
be easily diagonalized along with H\2. This is a first order 
degenerate perturbation theory. Because t± will turn out to be 
so small, we do not expect a large correction. 
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At each wave vector k in the RBZ, we need to diagonalize 
a 4 x 4 matrix to extract the energy eigenvalues. This matrix 
is 



( e k iW k t ± (k) \ 

-iW k e k+Q t ± (k+Q) 

tx(k) e k iW k 

t±(k + Q) -iW k e k+Q 



(5) 



The energy eigenvalues are, with s = ±1, 
2st x (k) ^ 



£k + £k+Q 



± ■ 



4W^ 



(6) 

From t± (k) it is clear that the electron pockets will be much 
more affected by it than the hole pockets. The bilayer splitting 
of the main frequency in the dHvA meausurement, presumed 
to be from the electron pockets, if it is to occur, should be 
smaller or of the order of the half width at the half maximum 
of the peak in the Fourier spectra, otherwise it would have 
been already resolved Q HOI . This fact combined with the 
Luttinger sum rule strongly constrains its magnitude. 

The parameters we choose for YBCO at 10% doping 
are 02: t = 0.3 eV, f = 0.3t, t" = f/9.0, t± = 8 meV, 
and Wq = 0.0825 eV. With these choices of the parame- 
ters and the chemical potential /i set to —0.2627 eV, we get 
the total hole doping of rih ~ 10%, constrained by the Lut- 
tinger sum rule |8), a sine qua non of the Fermi liquid picture. 
The corresponding dHvA frequencies according to the On- 
sager formula are F x « 944 T, F 2 « 967 T, F 3 w 570 T, and 
F 4 ps 450 T. The frequencies Fx and F 2 correspond to the 
hole pockets and are essentially the same within our accuracy, 
while F3 and F4 correspond to the electron pockets split by 
the bilayer coupling. The Fermi surfaces are shown in Fig.[T] 
Incommensurate order: In view of a recent experi- 
ment iflOl . it is interesting to pursue what incommensurate 
DDW order should predict for dHvA. Mean field theory with 
a class Hamiltonians show that it is difficult to achieve in- 
commensuration, unless the d-wave order parameter in the 
particle-hole channel is mixed with a s-wave component (may 
be natural in the presence of orthorhombicity) or a id xy com- 
ponent [16|. The reason is that it is energetically favorable 
to add holes at the nodes, unlike SDW, and provides a natural 
rigid band shift forming hole pockets and the consequent elec- 
tron pockets from the band folding. Nonetheless, it is possible 
that the Fermi surface could eventually move away from the 
nesting wave vector as a function of doping Q [TT). Whether 
or not this happens before the DDW gap collapses with in- 
creased doping is a very difficult question to answer. Any in- 
commensurate order is generally complex, especially if the in- 
commensuration is irrational (a slight abuse of terminology). 
Nonetheless, a useful approximation is to keep the largest gap 
and the hierarchy of gaps can be washed out due to thermal 
fluctuations, disorder, or magnetic breakdown fl8l . A simple 
approximation [ 1 3 1 is given by the particle hole condensate 
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FIG. 1: (Color online)Fermi surface split by bilayer coupling for 
commensurate DDW. The hole pockets centered at |(±7r/a, ±7r/a) 
are essentially unsplit by the bilayer coupling, In contrast, the elec- 
tron pockets centered at (n/a,0) and symmetry related points are 
split as described in the text. 



conserves current to only quadratic order in q, which may be 
sufficient for practical situations. 

There is some evidence from neutron scattering of incom- 
mensurate SDW fluctuations q = 7r(±2?7,0)/a and q = 
7r(0,±2?7)/a, where r] - 0.1 in underdoped YBCO 02- A 
similar estimate should also apply to DDW because the de- 
termining competition between the kinetic and the interaction 
energies are similar within a mean field theory; the precise 
value of j] is not particularly important at this time. We further 
assume a single wave vector q by spontaneous breaking of in- 
version symmetry, guided by the experimental observations. 
Note, however, that the product of inversion and time rever- 
sal, which is broken by DDW, is preserved, and we do not 
expect macroscopic currents. (One could have taken the or- 
der to have wave vectors q and — q, and thereby not breaking 
inversion symmetry, resulting in a different Fermi surface.) 
The excitation spectrum can be trivially solved. The Hamil- 
tonian is the same as in Eq. [5] except that the wave vector 
k now runs over the full Brillouin zone and W k must be re- 
placed by W k -> \{W k - W k -K_) and Q by Q -> K. The 
resulting Hamiltonian is easily diagonalized and the Fermi 
surfaces are shown in Fig. [2] The corresponding dHvA fre- 
quencies according to the Onsager formula are F± s» 1661 T, 
F 2 ~ 251 T, F 3 re 535 T, and F 4 ps 442 T. The frequen- 
cies Fx and F 2 corresponding to the hole pockets are essen- 
tially unchanged by the bilayer coupling, whereas F3 and F4 
corresponding to the electron pockets are split by the bilayer 
coupling. 

Robustness of dHvA in the mixed state: Many supercon- 
ductors exhibit dHvA effect deep in the mixed state (fields 
as low as half the upper critical field) with frequencies un- 
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FIG. 2: (Color online)Fermi surface split by bilayer coupling for in- 
commensurate DDW. The hole pockets centered at | (±7r / a, ±7r /a) 
are essentially unaffected by the bilayer coupling, The electron pock- 
ets are split as described in the text. The band structure parameters 
are the same as in Fig.[T] Wo = 0.115 eV, t± — 6 meV, r/ = 0.09, 
/i = —0.2575 eV, and nj, w 10%, as before. The chosen order pa- 
rameter spontaneously breaks time reversal symmetry. Thus, when 
reflection symmetry is also broken, an asymmetric band structure is 
not only allowed, but is expected on general grounds. 



shifted from the putative normal Fermi liquid state, except 
that the amplitude of the oscillations is diminished GUI . A 
particularly interesting analysis was given by Stephen ETI . 
To show that the same is true if the normal state is the 
DDW with a reconstructed Fermi surface, it is sufficient to 
examine the hole pockets close to the nodal points where 
the Hamiltonian is given by an effective Dirac Hamilto- 
nian, as the electron pockets are described by non-relativistic 
fermions and Stephen's results can be immediately taken over. 
The nodal Hamiltonian in the Landau gauge is ll22l Hp = 
^JvpVD [psO"3 + (Py — eBx/c)&2] , where cr's are the stan- 
dard Pauli matrices and the spatial anisotropy was removed 
with the redefinition of the coordinates x — \J vjj /vpx and 
V = \Zvf/vdD, where vp is the velocity perpendicular to a 
constant energy contour and vd is the velocity tangential to 
it. The momenta should similarly be rescaled, and we denote 
them by k. The unperturbed Green's function resulting from 
this Hamiltonian is used to solve the Gorkov equations l23l to 
find a self energy matrix, T, ni ^ 1 ,a 1 ;n 2 ,k2,a 1 (*^), which being 
diagonal in the Landau index and the wavevector (also inde- 
pendent of it), is, assuming that the Landau level index n ^ 1 
(the full expression can be found in Ref. [24 1), 
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i = \JhcfeB\ a = — 1 for the lower branch of the Dirac 
spectra and a = +1 for the upper branch . For simplicity we 
have used an s-wave superconducting gap A, but it is straight- 
forward to generalize to a d-wave gap. The main result is un- 
changed, namely that the real part is very small compared to 
the imaginary part. (An interesting aspect of this formula is 
that only for the component ct\ = a 2 = — 1 the result is non- 
vanishing.) Thus, we can define an "impurity" scattering rate 
at the Fermi energy due to vortices in the mixed state as 



A 2 



(8) 



(7) 



where the effective mass for the DDW quasiparticles is de- 
fined by m* = /i 2 / (vpVjj) and lj* = eB/m*c. In agreement 
with Stephen, the cyclotron motion is so much faster than the 
vortices that the electrons see the vortices as static impurities. 
It is now obvious that there would be an additional Dingle fac- 
tor arising from the vortices in the mixed state, but the dHvA 
frequencies will be unshifted to an excellent approximation 
from the putative normal DDW state; the superconducting gap 
sets only the magnitude of t v . A more complete discussion 
will be given in Ref. [24|. It is not difficult to see, however, 
that for a d-wave superconducting gap A 2 will be replaced by 
the Fermi surface average of the square of the d-wave super- 
conducting gap. 

Coda: There are important questions that remain unan- 
swered. The most glaring is the Fermi surface determined in 
angle resolved photoemission spectroscopy (ARPES), which 
is mostly consistent with Fermi arcs [25] in the relevant dop- 
ing regime in the normal state. It is noteworthy, however, that 
in electron doped cuprates, the evidence of hole and electron 
pockets is a solid but an often forgotten fact \26\. More re- 
cent ARPES experiments have indicated remarkable data for 
hole pockets 1271 . In contrast, experiments in YBCO films 
where the doping is adjusted by applying potassium overlayer 
to reduce it to 10% see once again a Fermi arc [28 1. Although 
by itself this is not in contradiction with the notion of a hole 
pocket because it has long been argued that the back side of 
the hole pocket should have very small intensity in ARPES 
simply from the DDW coherence factors [29|. The real issue 
is of course why the electron pockets are not seen, at least the 
two sides of it. Perhaps adding potassium overlayer may cause 
significant charge scattering and the electron pockets may be 
highly fragile with respect to disorder. 

One might argue that the ortho-II potential of YBC06.5 
might be important [30|, a material for which data for Hall, 
SdH, and dHvA data exist. We believe that this may not be 
so because highly polarizable BaO-layers next to the chains 
should screen the potentials quite effectively, to the extent that 
even disordered chains appear to have little effect in the pla- 
nar physics in many properties. It is also important to rec- 
ognize that SdH (and dHvA) measurements are also available 
in YBa2Cu4C>8, a double chain compound, and even nega- 
tive Hall coefficient is clearly observed in agreement with the 
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ortho-II materials. It is difficult to believe that this universality 
is achievable if the chain potentials were playing a significant 
role. Note, however, that as yet no oscillations are observed 
in ortho-VIII, although negative Hall coefficient is similarly 
observed Q. 

We are fully cognizant that neutron scattering experiments 
have resulted in a controversy, as commented in Ref IU~2l . with 
regard to DDW. The fact remains that if the pseudogap is sup- 
posed to reflect an order with a strongly momentum depen- 
dent gap in the particle-hole channel, this is impossible with- 
out breaking translational symmetry. 

Bilayer splitting in overdoped materials is well estab- 
lished [31 1, although the actual splitting, 88 meV, is dramat- 
ically smaller than the local density band structure calcula- 
tions, which is 300 meV. The reason for this is likely to be 
the Z-factor associated with the Green function |23l . We have 
argued that the splitting at doping of 10% is even smaller and 
only about 12-16 meV. If this is the case, a more precise un- 
derstanding of this renormalization is necessary, but it could 
not be a breakdown of the Fermi liquid theory, as these oscilla- 
tion measurements are indicative of quasiparticles of a Fermi 
liquid. 

So far the high frequency peak at 1654 T observed in DC 
measurements ifTUl has not been observed in pulsed field mea- 
surements 0. The commensurate order, as we have seen, pre- 
dicts three main frequencies, two of which are from bilayer 
split electron pockets, but the hole pocket remains essentially 
unsplit at about 970 T; however, this hole pocket frequency is 
not yet seen (although there is a very weak signature in some 
data), may be because it lies too close to the second harmonic 
of the electron pocket frequency. On the other hand, if the or- 
der is incommensurate, the observed frequency at about 530 
T (electron pocket) should be similarly split by the bilayer 
coupling. The hole pocket frequency at about 1650 T will es- 
sentially remain unsplit because of the nature of the bilayer 
coupling. To conserve the Luttinger sum rule, however, there 
should be a lower hole frequency at about 250 T, whose split- 
ting will be similarly very small. 

Although there are other analyses of the quantum oscil- 
lation phenomenon l32l l33l l34l [35], the prediction of the 
splitting of the main frequency due to bilayer coupling and 
a smaller hole pocket frequency at 250 T for the incommen- 
surate case are entirely new, and, if verified, would constitute 
a striking argument for the Fermi liquid picture — for a doping, 
as low as 10%, this would be surprising to say the least. 
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Note added: After our work was completed, we received 
a complementary paper ll36l . It emphasizes ortho-II ordering 
that results in a very different Fermi surface. We have ar- 
gued above that ortho-II potential is likely to be unimportant. 
The treatment of bilayer coupling in Ref. [36| does not yield 
the splitting that we have emphasized. Moreover, neither the 
incommensuration nor the analysis of the mixed state is dis- 



cussed there. 
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